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Abstract 

Analytical expression for energy of eigen-modes in magnetohydrodynamic flows of ideal fluids 
is obtained. It is shown that the energy of unstable modes is zero, while the energy of stable 
oscillatory modes (waves) can assume both positive and negative values. Negative energy waves 
always correspond to non-symmetric eigen-modes - modes that have a component of wave-vector 
along the equilibrium velocity. These results suggest that all non-symmetric instabilities in ideal 
MHD systems with flows are associated with coupling of positive and negative energy waves. As 
an example the energy of eigen-modes is calculated for incompressible conducting fluid rotating in 
axial magnetic field. 



1 



Energy consideration is of primary significance in stability analysis of different magneto- 
hydrodynamic (MHD) systems. It is well known that the energy associated with the waves 
(purely oscillatory eigen-modes) may change its sign and become negative . The en- 
ergy should be withdrawn from the system to let the negative energy wave be excited. So, 
a negative energy wave is a potential source of instability since no extra energy is needed to 
increase its intensity. Instability can arise, for example, if a negative energy wave is subject 
to external dissipation; then the subsequent removal of energy from the wave will cause it to 
grow. In a conservative system, the instability can occur due to the simultaneous excitation 
of positive and negative energy waves. In this case, energy is transferred from the negative 
energy wave to the positive energy wave, allowing both modes to grow and the total energy 
to remain constant. Waves having energies of various signs enable researches to explain 
different types of instabilities in fluid dynamics j^. 

In the present paper we calculate the energy of the eigen-modes in ideal one-fluid MHD 
and show that all instabilities of non-symmetric eigen-modes in MHD systems with equilib- 
rium flow are related to the coupling of negative and positive energy waves. Following Ref. 
J], we consider linearized dynamics of displacement vector ^ 

P|| + 2P(V-V)f -F[C]=0, (1) 

where p and V are stationary values of fluid density and velocity, respectively. The general 
form of linearized force operator F[^] in ideal compressible MHD is 

m = -p(V-V)2^ + p(^-V)(V-V)V + 

+ V- (p^)(V- V)V- V(5P+ (2) 

+ ^ (V X 5B) X B + — (V X B) X 5B. 

An An 

Here, B is equilibrium magnetic field and 

(5B = V X (I X B) 

is its perturbation. The perturbation of fluid pressure SP can be specified by thermodynamic 
properties of the system. For example, if the process is adiabatic with adiabatic index 7 
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then 

5P = -I . VP - 7PV • ^. 

In the case of incompressible MHD, such equation appears to be excessive, instead one has 
to impose the incompressibihty condition V ■ ^ = 0. 

A number of formal properties of Eq. ([1]) can be established. Force operator F[^] is 
Hermitian (self-adjoint) in the following sense, 

j ri ■ F[|] dh = j F[t7] dh, (3) 

while the second term in Eq. ([1]) is antisymmetric: 

j prj-{V-\/)^dh = - j p^- {V ■V)ri£r. (4) 

Integration in Eqs. ([3]) and (jlj) is performed over the fluid volume under the assumption 
that displacements on the boundary vanish. 

In our subsequent discussion, the displacement vector ^ is supposed to be complex. In 
order to obtain the correct expression for energy of perturbations in this case, we multiply 
Eq. ([T]) by complex conjugate d^*/dt and integrate over the space: 

The complex conjugate of this equality is 

Summing up these two equations and using the properties ([3]), (j4]) we arrive at the energy 
conservation law in the form dE/dt = 0, where 



dt 



e ■ m dh. (5) 



As usual in mechanics, the total energy of the perturbations consists of kinetic part (first 
term) and of potential part (second term). 

Since the equilibrium quantities have no time dependence, we look for a normal-mode 
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solutions to Eq. ([T]) in the form 

|(r,t)=^(r)e-'^*. (6) 
Then, the equation of motion ([1]) leads to eigen-value problem 

uj'^pi + 2iujp{\ ■ V)i + F[^] = 0. (7) 

Multiplying this equation by complex conjugate ^ and integrating over the fluid volume, 
we arrive at quadratic equation for eigen- frequency a;, 

Auj'^ -2Buj-C = (8) 

with coefficients 



A = j pl^prf^r > 0, 

B = -I j pi ■ (V- V)^rf^r, 

C = - J k*F[^]dh. 



Solving Eq. ([8]) we find 



B ± VB^ + AC 
^1,2 = J ■ (9) 



This expression allows to determine eigen-frequency corresponding to known eigen-mode ^. 
Since all coefficients in Eq. ([8]) are real [for coefficients C and B it follows immediately from 
properties and (jlj), respectively], the instability in the system is possible if and only if 
B^ + AC < for some eigen-mode. 

Now we are able to determine the energy of the eigen-mode with eigen-frequency IQ. 
Substituting ([6]) into expression ([5]) we obtain: 

E=^-{A\u;\' + C). (10) 



In the case of unstable mode, B^ + AC < 0, so 



TABLE I: Eigen- frequencies oo 1,2 and corresponding energies £'1,2 of stable eigen-modes for different 
values of coefficient C (i? > is assumed for simplicity). 
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and the energy is 



Ei,2 = 0. 



For stable mode, B^ + AC > and the energy is 



i?i,2 = ^^'^"^^ {^/WTAC±B 



±ujVW+AC. 



(11) 



(12) 
(13) 



Therefore, energy of stable eigen-mode changes the sign if its frequency changes the sign. 

Depending on the system parameters different options are realized in the case of stable 
eigen-modes (Table [I]). As one can see, there is an interval of parameters at which the 
waves with positive and negative energy coexist (option 2). One boundary of this interval 
corresponds to the stability threshold (option 1), the other - to change of sign of eigen- 
frequency UJ2 (option 3). This result suggests that the instability in the ideal MHD system 
with flow can be associated with coupling of positive and negative energy waves. 

We note here that all negative energy waves are non-symmetric modes, i.e., they have 
spatial dependence along the equilibrium flow, so, the coefficient B ^ 0. For symmetric 
modes or in the absence of flow we have B = and the energy is 



E = -i\C\+C). 



Therefore, energy of symmetric modes is never negative, and their stability can be inves- 
tigated by use of energy principle jj]. In a case of non-axisymmetric modes, the energy 
principle fails and special arrangements should be made to modify it (see, e.g., j^). 

In order to verify the above analytical results, we calculate the energy of eigen-modes 



5 



of incompressible fluid rotating in homogenous transverse magnetic field B = Be^. The 
equilibrium velocity profile used in our calculations corresponds to the electrically driven 
fiow in circular channel and has a form 

V = rfi(r)e^, ^](r) = ^ (14) 

in cylindrical system of coordinates {r, (f, z}. Here, ri and r2 are inner and outer radii of the 
channel, respectively, and Qi is the angular velocity at ri. This type of fiow profile is used 
in new experimental device for laboratory testing of the so-called magnetorotational 
instability (MRI), which plays an important role in many astrophysical applications (see 
reviews |7|, l8| ) . 

A detailed eigen-mode analysis of such fiow has been performed in Ref. Assuming 

one obtains eigen- value problem 

i^ + ™^^ + ,t.j.=0 ' (15) 

with boundary conditions 

where 

is Alfven frequency. 



Un) = Ur2) = 0, (16) 



KB 



uo = uo — rriQ 

is "shifted" eigen- frequency and 611 is perturbation of the total normalized pressure, 

P 

n = - + — . 

p Airp 



6 



A general expression for energy of perturbations ([5]) for this system reads 



E=P- 
2 



dt 



or 



(17) 



Substituting ^ from the eigen- value problem (fTSjl . (|T6|l we find the energy of stable mode 
with frequency u: 



E 



nphoj I rluj 



+ 



(18) 



ri 



+ 



+ 



where h is the height of the channel. For axisymmetric eigen-modes with m = this 
expression is reduced to 



r2 



E = Tiphoj'^ 



klr"^ 



djr^r 
dr 



+ 



r dr. 



Therefore, the energy of axisymmetric eigen-modes is always positive if 7^ 0. Formally, 
this case is described by f lT2|) with coefficient B = 0. 

In Figs. [H [2] the calculated dependencies of frequency and energy for two potentially 
unstable eigen-modes on the parameter Qi/uja are shown. In the axisymmetric case (m = 0), 
both branches of energy are positive and coincident (Fig. Wp)- The merging point in Fig. 
[It corresponds to Qi/uja ~ 2.0 which is the threshold of magnetorotational instability for 
m = 0. The nature of axisymmetric MRI is not related to the subject of negative energy 
waves and can be explained by the mechanism similar to one of Raleigh- Taylor instability 

Q. 

For m = 1 the behavior of both energy curves in Fig. [2)d is completely described by Eq. 
ffT2]) . The MRI threshold in this case is Qi/uja ~ 1-7. When 1.1 < Qi/ua ^ 1-7 the positive 
and negative energy waves can coexist in the system. At Qi/uja ~ 1.1 the frequency uj2 
changes the sign (Fig. [2^, dashed line), so both energy branches become positive. 

It should be noted that merging points in Figs. [1] and [2] determine the magnetorotational 
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FIG. 1: Calculated dependence of eigen-frequency (a) and energy (b) on ratio ^i/lva for two most 
unstable eigen-modes with azimuthal number m = 0. Energy is given in arbitrary units. 



instability threshold. In the flow given by fll4p this threshold decreases with azimuthal 
number m, as discussed in Ref. jsl. For large m it approaches the asymptote 



VLi _ 2 

Ua m{l-ryrl)' 



(19) 



The calculated dependence of MRI threshold on small m is presented in Fig. [3l 

In conclusion, we have shown that all non-symmetric MHD instabilities in ideal fluids 
with flows can be explained as a coupling of originally stable positive and negative energy 
waves. These results are supported by calculations of frequencies and energies of eigen-modes 
in the flow that can be unstable with respect to magnetorotational instability. 
This work is supported in part by NSERC Canada. 
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FIG. 2: Calculated dependence of eigen-frequency (a) and energy (b) on ratio Qi/oja for two most 
unstable eigen-modes with azimuthal number m = 1. Energy is given in arbitrary units. 




FIG. 3: Calculated dependence of magnetorotational instability threshold on azimuthal mode 
number m (solid line) and its asymptote for large m (dashed line). 
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